Heat capacity and phonon mean free path of wurtzite GaN 
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^S) . We report on lattice specific heat of bulk hexagonal GaN measured by the heat flow method 

in the temperature range 20-300 K and by the adiabatic method in the range 5-70 K. We fit the 
Oj' experimental data using two temperatures model. The best fit with the accuracy of 3 % was obtained 

.^fH , for the temperature independent Debye's temperature Ou = 365 K and Einstein's temperature 

^ ■ 6e = 880 K. We relate these temperatures to the function of density of states. Using our results for 

CN) ' heat conduction coefficient, we established in temperature range 10-100 K the explicit dependence 

of the phonon mean free path on temperature /ph oc T"^. Above 100 K, there is the evidence of 
contribution of the Umklapp processes which limit phonon free path at high temepratures. For 
phonons with energy fee x 300 K the mean free path is of the order 100 nm. 
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^ , PACS numbers: 65.40.-|-g;63.20.Dj 
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I. INTRODUCTION 



High pressure grown GaN has been used as a substrate for high quality epitaxial layers, heterostructures and blue 

laser diode J,] . Efficient heat removal is critical to the performance of high power semiconductor electronic devices, 

■ O ■ such as light emitting diodes, lasers, photodetectors or metal-semiconductor field effect transistors, and high electron 

mobility transistors. For this reason, thermal conductivity and heat capacity are mandatory characteristics of GaN, a 

very topical semiconductor these days. In our previous paper we studied thermal conductivity of this compound M- 

Early measurements of the heat capacity of GaN had been performed for polycrystalline powder samples Q]- 

Measurements of the heat capacity of hexagonal single crystals of GaN in the 20-1400 K temperature range were 

\j| ' reported by Kremer et al Q and for 320-570 K range for powder samples by Leitner et al [a|- Their results match 

J!Z. findings of Kashchenko et al for powder samples [31 very well. 

r — In our paper we focus our attention on the bulk wurzite GaN heat capacity Cp measured in the low temperature 

fSl ! region. Having measured Cp and thermal conductivity coefficient k |2|, we are able to estimate the temperature 
^^ ■ dependence of the mean free path Zph of acoustic phonons. The knowledge of Zph is important for understanding the 
^D ■ thermalization processes in GaN substrates and structures as well as for the heat thermal budget calculations. 
>0 , The phonon density of states (DOS) g{iLj), which includes contributions from all phonons over the entire Brillouin 

zone, is needed for the calculations of various thermodynamic characteristics, e.g. heat capacity and thermal con- 
ductivity. As neutrons probe phonon modes throughout the Brillouin zone, the neutron scattering can be applied 
to establish DOS. Nipko et al experimentally estabhshed DOS for powder samples of GaN [g. The partial DOS of 
S I near zone center optical phonons can be also obtained using Raman and infrared spectroscopy (in the case of GaN 
r^ . cf. references given in Jg|). DOS of GaN was calculated by Nipko et al p, Azuhata et al 0, Bungaro et al [g and 
P5 ■ Sanati and Estreicher [3. 

O ■ Debye proposed a simple form of DOS based on the elastic continuum theory 

O 

JH , where 9{x) is the Heaviside step function, gui'-^) depends only on one parameter - the Debye frequency o^q (or on 
the Debye temperature 0d = fii^D/kB)- The Einstein form of DOS 

5E = 3NsSiuj - we) (2) 

where S{x) is Dirac's delta distribution, N - the number of elementary cells, and s - the number of particles in it, 
assumes that all lattice particles vibrate with the same frequency of optic phonons, lue, hence the DOS depends also 



^2 

.gD(w) = 3^r-0{ujB - LO), (1) 
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FIG. 1: Dependence of heat capacity Cp'^ on temperature measured by two methods: adiabatic calorimetry in the temperature 
range 5-70 K and the heat-flow method in the range 20-300 K. Thin line represents results for the dual-temperature model. 
Open circles and squares represent results of measurements by Kremer et al p|. Thick line - results of our measurements. The 
contributions of Einstein's (curve no 1) and Debye's models (curve 2) are also shown. 



on the Einstein frequency lu-e or on the corresponding temperature 9^ = fiw-pJk-Q. Debye's and Einstein's models with 
several characteristic temperatures (frequencies) were also introduced (cf. [lOJ). 

In our paper we discuss the results of measurements of the wurzite GaN heat capacity in terms of a model depending 
on the Debye and Einstein temperatures [13 , namely we assume that Cp [ J/mol K] 



Cp(T)«CD(r) + CE(r) = 

37VAA:B3(r/0D)'F(T/0D) 



(3) 



+37VAfcB [e^m 



2 -Se/T / (gSE/T 



where 



FiT/Oj,) = 



T/en 



dxx'^eyie'' -if 



(4) 



and -/Va is the Avogadro number. 

Hexagonal GaN single crystals were grown at the Institute of High Pressure Physics, Polish Academy of Sciences 
[11] and in ETH Ziirich. The specific heat of approximately 30 mg GaN single crystal was measured using the heat- 
flow calorimeter |l2| . In this method, Cp measurements could be performed both upon cooling and heating. The 
high data density is another advantage of this method (approximately 10 points/K). Using this method we measured 
heat capacity in the temperature range 20-300 K. Adiabatic heat capacity measurements were performed in the 
temperature range 5-70 K on three different samples with masses approximately 100 mg. 

In Fig. n full circles represents the total specific heat of GaN single crystal measured by heat flow method upon 
heating and cooling in the temperature range 20-300 K and by adiabatic method in the range 5-70 K (thick solid line). 
In the common temperature interval 20-70 K, the results obtained by both methods coincide. 

We consider the contributions of the Debye heat capacity Cd and the Einstein heat capacity Ce to the total heat 
capacity Cp ~ Cd + Ce. For temperature independent Debye's and Einstein's temperatures 6b — 365 K, 9e — 



K, the plot of the function (0J approximates the heat capacity obtained in our experiments with 3% accuracy (Fig. 
^. The Cp values measured below T = 7 K are less accurate. If Cp exhibits the T^ behavior, and if the free carriers 
contribution is negligible, then a Cp/T versus T^ plot results in a straight line through the origin. This means that 
for the specimens studied, the electron contribution to Cp is negligible. Fig. |21 shows Cp/T versus T^ plot in the 
range of 0-45 K. From this figure it can be concluded that Cp/T is proportional to T^. 

We also calculated the Debye temperature using familiar formulas 13] and the room temperature elastic constants 
for GaN. The results are summarized in TableQ] The value of Ou providing the best fit of Cp(T) is smaller than values 
collected in Table ^ It is interesting to note that in the case of measurements for powder samples J3| O^ depends on 
temperature. 

The main feature of the phonon DOS of hexagonal GaN [a|-[3 is the presence of two well separated bands of low 
energy phonons (acoustic and optic) in the range 0-42 meV and optic phonons in the range 70-90 meV (cf. Fig. [SJ. 



o 



0.15 



0.10 



0.05 



0.00 




r 



-t(cxp) 



FIG. 2: Cp /T versus T for hexagonal GaN in the low temperature region. The sohd line presents Debye's heat capacity 
Cd with 9u ~ 365K. The experimental points are shown as open circles. The squares represent experimental results obtained 
by Koshchenko et al Q for powder GaN samples. 
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TABLE I: Debye's temperature for various sets of elastic constant and our result 

^) C. Degcr, E. Born, H. Angerer, O. Ambachcr, M. Stutzmann, J. Hornsteinor, E. Riha, G. Fischcraucr, Appl. Phys. Lett. 72, 2400 (1998) 

^) V. A. Savastcnko, A. U. Shclog, Phys. Stat. Sol. A 48, K135 (1978) 

^) A. Polian , M. Grimsditch , I. Grzegory, J. Appl. Phys. 79, 3343 (1996) 

*) Y. Takagi, M. Ahart, T. Azuhata, T. Sota, K. Suzuki, S. Nakamura, Physica B 219/220, 547 (1996) 

^) M. Yamaguchi, T. Yagi, T. Azuhata, T. Sota, K. Suzuki, S. Chichibu, S. Nakamura, J. Phys. C 9, 241 (1997) 

'^) R. B. Schwarz, K. Khachaturyan, E. R. Weber, Appl. Phys. Lett. 70, 1122 (1997) 

'') T. Dcguchi, D. Ichiryu, K. Toshikawa, K. Sckiguchi, T. Sota, R. Matsuo, T. Azuhata, T. Yagi, S. Chichibu, S. Nakamura, J. Appl. Phys. 86, 

1860 (1999) 

**) M. Yamaguchi, T. Yagi, T. Sota, T. Dcguchi, K. Shimada, S. Nakamura, J. Appl. Phys. 85 , 8502 (1999) 

Therefore, analysis of our Cp(T) data with the help of two temperatures 9^ and 8e is sound. The cut-off energy 
^B^D = 32 meV is located near the upper Hmit of the lower energy band. The Einstein energy /cb^e = 77 meV is 
located near the center of the energy band for optic phonons. 

Consider the kinetic expression for the thermal conductivity coefficient k 



K = Cplphv/3, 



(5) 



where Iph is the phonon mean free path and v is the mean velocity of acoustic phonons. Using it we estimate the 
experimental phonon mean free path Zph- Obviously v and Iph are characteristics of the "propagative" phonons, i.e., 
those which really contribute to the transport phenomena. Consequently the heat capacity Cp to be considered should 
be only due to the contribution of these phonons, i.e., Cp « Cd- For GaN 5 « 5 x lO^cm/s. Since we have established 




FIG. 3: Ab initio phonon dispersion curves and density of states DOS for GaN in the wurtzite structure 
is the Debye DOS Eq. Q for huio ~ 32 meV. The Einstein energy huiE = 77 meV is also indicated. 



The dashed line 



the temperature dependence of k U of the specimen which was also used in our heat capacity measurements, we can 
find the dependence of Zph on temperature. The results of such analysis are presented in Fig. |21 Below 100 K the 
mean free path is proportional to T^^. Above 100 K the Umklapp-processes gain importance 14], making mean free 
path shorter. For 300 K the phonon mean free path exceeds 100 nm, which is 10 times greater than characteristic 
length of nanodcviccs based on low dimensional gases of carriers |15j. 16J. This means that in such devices the energy 
removal by acoustic phonons with energies lower than fce x 300 K is ballistic. This observation is in agreement with 
our previous results for heat pulses propagating in bulk GaN \u\ . 
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FIG. 4: The temperature dependence of the phonon mean free path for the wurzite GaN monocrystalline specimen. Triangles 
(line no. 



/,-ir'(<=''p) 



1) indicate V-"^^^' = 3Kcxp(r)/'uCp ^ (T) and squares (line no. 2) stand for /ph = 3Kexp(r)/iiCD(T). The broken line 



Ph 



represents results obtained by Liu and Balandin [Tj|. Insert - heat conductivity coefficient IJ- 
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